Fast generation of spin-squeezed states in bosonic Josephson junctions 
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We describe methods for fast production of highly coherent-spin-squeezed many-body states in 
bosonic Josephson junctions (B J Js) . We start from the known mapping of the two-site Bose-Hubbard 
(BH) Hamihonian to that of a single effective particle evolving according to a Schrodinger-like 
equation in Fock space. Since, for repulsive interactions, the effective potential in Fock space 
is nearly parabolic, we extend recently derived protocols for shortcuts to adiabatic evolution in 
harmonic potentials to the many-body BH Hamiltonian. The best scaling of the squeezing parameter 
for large number of atoms A'' is ^| ~ 1/N. 
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Recent experimental results from the group of M. 
Oberthaler have provided evidence of the generation of 

fin-squeezed states in bosonic Josephson junctions [ij- 
Two different setups with ultracold ^''Rb bosons 
have been explored: a) a cloud trapped in an exter- 
nal double-well potential [l], Q (external BJJ) and b) a 
cloud of atoms in two different hyperfine states trapped 
in a single harmonic potential with a linear coupling 
between the two internal states [i] (internal BJJ). In 
both cases, the atomic-clouds are well described by the 
two-site Bose-Hubbard model |5|-|7[, a Lipkin-Meshkov- 
Glick-like (LMG) Hamiltonian [8|. As previously pointed 
out 0, this Hamiltonian can produce largely squeezed 
spin states and, correspondingly, large sets of entangled 
pseudo-spins d, [13, U | • 



the systerr0 on the number of atoms, N, the atom-atom 
interactions, U, and the linear coupling J. To do so, we 
benefit from the explicit mapping, in the large N limit, 
between the two-site Bose-Hubbard model and an effec- 
tive single particle system described by a Schrodinger-like 
equation |15l4l7| (in the more general context of the LMG 
model, this connection was established earlier 3 3: al- 
beit with no known relation at the time with ultracold 
atom physics). 

The dynamics of a bosonic Josephson junction can be 
well described by a quantized two-mode model [i-Q. For 
the external BJJ case 0, , or by construction in the 
internal BJJ case, the system of ultra-cold bosons is well 
modeled by the Bose-Hubbard Hamiltonian UHbh 



n 



BH 



~2JJ^ + UJi 



(1) 



Simultaneously, and due to the need to control and im- 
plement quantum resources, there has been an increasing 
interest in developing fast protocols to shortcut adiabatic 
following [l2l . [l3| . The key purpose is to engineer pro- 
cedures to drive, in a finite time, a system from some 
initial state to a final state that could be reached with 
an adiabatic, slow process. The different protocols may 
in addition aim at a) minimizing the transient energy 
excitations, b) reducing the sensitivity to noise, or c) 
minimizing other variables of interest 1^. It should be 
noted that, in general, the initial and final states arc not 
necessarily required to be the ground states (GSs) of the 
system. In this work, however, we concentrate in tran- 
sitions between ground states corresponding to different 
values of the model parameters which can be controlled 
externally. The proposed method is designed so that the 
desired final state is produced as a stationary eigenstate 
of the Hamiltonian, with no need to stop or freeze the 
dynamics. Analytical formulas to perform this type of 
processes exist for the harmonic oscillator [l2| . We show 
that they can be extended to many-body states described 
by the BH Hamiltonian. That allows us to compute the 
dependence of the attainable coherent spin squeezing in 



where the pseudo-angular momentum operator J = 
{Jx, Jy, Jz\ is defined as, Jx = (l/2)(a|a2 -I- 6201), 
Jy = {l/2i)(a\a2 — a\ai), and Jr = (l/2)(a|ai — 020.2)- 
dj creates a boson in site j, [di,dj] = 5i^j. J is the hop- 
ping strength, taken positive, and U is the non-linear 
coupling strengthproportional to the atom-atom s-wave 
scattering lengtfo. In this work we consider repulsive 
interactions, C/ > 0. For internal BJJ, the inter-species 
s-wave scattering length in ®^Rb atoms can be varied 
by applying an external magnetic field thanks to a well 
characterized Feshbach resonance &t B = 9.1 G, as done 
in Ref. [3| for the setup that we are considering, thus 
permitting precise control over the U . In this work, we 
use a time dependent U{t), keeping J and N fixed dur- 
ing the time evolution. To do so one needs to track the 



^ The use of the spin squeezing as a witness of A'^— particle entan- 
glement is established in Ref. [llj . In this letter we focus only in 
the characterization of the squeezing. 

^ In the internal BJJ, U is proportional to ai,i +02,2 — 2ai_2, with 
ai^i and 02,2 the intra-species scattering lengths and ai^2 the 
inter-species one 
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resonance of the linear coupling (responsible for the hop- 
ping term in the Hamiltonian) which is altered due to 
the second order Zeeman shift when B is varied. Albeit 
technically challenging, this is within reach of the current 
experimental setups [20[. 

The time dependent Schrodinger equation (TDSE) is 
written as, tdt\'^) = HbhI^')- For a given N, an ap- 
propriate many-body basis for this bosonic system is 
the Fock basis, {|m^ = (iVi - N2)/2)}, with = 
—N/2,...,N/2. A general many-body state, |^), can 
then be written as = J2mX-N/2 ^^n, \mz) . 

As customary, the number squeezing parameter is de- 



fined as @,[2l|2^, £_l{t) 



AJ^/(AJ^)bin where AJ^ = 



(Jf) — {Jz)"^ and (AJj')bin = N/i in the binomial case. 
The many-body state is said to be number-squeezed if 
< 1 [9[. The coherent spin-squeezing parameter is 
defined as i E E| = N{AJ^)/{X)^ = 
where the phase coherence of the many-body state is 
a{t) = {^{t)\2X/N\^{t)) . takes into account the 
delicate compromise between improvements in number- 
squeezing and loss of coherence. States with < 1 have 
been proposed to be used in a new Ramsey type atom in- 
terferometer with an increased phase precision compared 
to the coherent spin state. This gain in precision can be 
directly related to entanglement in the system [loj . 

As detailed in Refs. |15l4l7| a systematic expansion of 
Hn = 'H-bh/{NJ) in powers oi h = 1/N gives in the 
semiclassical limit, iV >> 1, the pscudo-Schrodinger-likc 
equation. 




FIG. 1. (color online) Time evolution of the coherence a{t) 
and the number squeezing of the system computed using 
the fast adiabatic evolution of Eq. Q (solid) compared to 
the case of linear ramping of A{t) (dashed). The correspond- 
ing instantaneous adiabatic results for both the fast adiabatic 
case and a linear ramping are given in dotted and dot-dashed, 
respectively. A(0) = 0, A{tB.M) = 20, and = 100. 



Fast shortcut to adiabatic evolution 



nNiz)i^{z) = -2h^d,Vl - z^d,^{z) + V(z)V'(^) (2) 



where V(z) = -Vl - + {1/2) Az^, z = m^/{N/2) , 
A = NU/{2J) and 4'{z) ~ ^jN/2cm^, normahzed as 
j^^dz\%l){z)\'^ = 1. The corresponding TDSE writes as 
ihdtip{z,t) = HAr'0(z,i) with time measured in units of 
1/J. In these units the Rabi time is defined as tRabi ~ 
tt/J. 



As ex plai ned in jl7l |. and previously noted by other 
authors (l5L [23l | , for repulsive interactions the potential 
in Fock space, V(z), is to a very good approximation a 
harmonic oscillator: Neglecting the z dependence of the 

zV2 



effective mass term and expanding vl 
in V{z): Eq. ^ reduces to 



1 



n 



N 



-2h'di 



1 2 2 
-OJ z . 



(3) 



with Lu = 2\J\ -I- A. 

We have checked that all the expectation values of Ji 
and jf, i = x,y,z computed by directly solving the 
TDSE for the BH Hamiltonian, and the corresponding 
ones using the solution Eq. ([2]) and the explicit expres- 
sions given in Ref. [l3| agree almost perfectly in all cal- 
culations reported in this letter. The only minor discrep- 
ancy is explicitly shown in Fig. [3] 



Our method for fast adiabatic-like preparation of a 
given ground state requires control of the atom-atom in- 
teraction at time scales of fractions of the Rabi time. In 
the parabolic approximation discussed above the proto- 
col developed in Ref. [1^ has a direct translation into our 
variables. We need to impose the following time varia- 
tion of aH, A{t) ^ w^(t)/4 - 1, where a;^(t) obeys the 
Ermakov equation, 



b{t)/J^ +LJ^{t)b{t)^Oj'^/b^t) 



(4) 



with time in units of 1/J as explained above. The poly- 
nomial ansatz for b{t) from Ref. [l2j takes the explicit an- 
alytic form 6(i) = 6{r -l)s^- 15{r-l)s^ + 10{r-l)s^ + l, 
s = t/tf, and r = ^Ju)Q|wf where wq = 2-^/1 -|- Aq and 

= VI + A/ • 

Our calculations confirm that this procedure, although 
not exact, is extremely accurate in the many-body simu- 
lations despite the fact that it was derived for the single 
particle problem in a harmonic oscillator Hamiltonian. 
This is due to the validity of the explained mapping of 



^ Note that the derivations in the previous section imply a fixed J 
and N , thus the variation in A has to be due to a variation in U . 
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FIG. 2. (color online) Considered A{t) (upper panel): three 
final times are considered, tp = 0.2,1 and 5tRabi. Fidelity 
(middle panel) when tp = 0.2 tnabi, and coherent spin squeez- 
ing ^1 (lower panel) as a function of t using the fast-adiabatic 
process described in the text. The adiabatic ^1, correspond- 
ing to the instantaneous ground state for each A, is shown in 
dashed lines in the lower panel. The initial and final A are 
and 200. N = 100 



the BH to the single particle pseudo-Schrodinger equa- 
tion, which turns out to be well approximated by a har- 
monic oscillator potential. This harmonic approximation 
improves as N is increased because the spread of the GS 
wave function ^ 1 /y/N, thus exploring only the very cen- 
tral part of V(z) in 2: e [—1, 1]. Exceptions appear in the 
ultra-fast preparation as discussed later. 

Let's start by comparing the fast protocol to shortcut 
(SC) the adiabatic following presented above to the case 
of a linear ramping A,(t) = A(0) + (t/tF)mtF) - A(0)). 
In this case we consider a relatively long final time of 
one fuU Rabi period iRabi = t^/J, see Fig. [T] For both 
A(<) we present the results of the corresponding instanta- 
neous ground state of the many-body system for compar- 
ison. First, we note that the SC does produce a final re- 
sult which has similar coherence a and coherent number 
squeezing as the corresponding adiabatically evolved 
case. In contrast, the linearly ramped system deviates 
notably from the adiabatically evolved state in the final 
time. As we arc considering fairly long times, the SC is 
not found to deviate appreciably from the instantaneous 
adiabatic following of the state, except for a short pe- 
riod at early times. The transient non-adiabaticity will 
be shown to increase as we require shorter tp^s. 

The fidelity between the evolved many-body state 
\^{t)) and the corresponding instantaneous ground-state, 
|(\E'(t)|vI'Gs)|, is plotted in Fig [2] and is seen to be ex- 
tremely close to 1 at the final time for all considered 
tp =0.2, 1 and 5 iRabi- In the figure we compare the 
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FIG. 3. (color online) Considered A{t) (upper panel), fidelity 
(middle panel) and coherent spin squeezing ^5 (lower panel) 
as a function of t for the exact Bose-Hubbard calculation 
(solid) and for the Schrodinger-like equation The final 
time tp is 0.05 tRabi. The initial and final A are and 500. 
N = 100. The inset zooms on the times close to tp. 



full many-body evolution to the corresponding instanta- 
neous ground state for different tp. As mentioned above, 
for short times, the fast-adiabatic passage produces inter- 
mediate many-body states which depart from the instan- 
taneous ones, as seen clearly in the drop of the fidelity 
for very short times in the tp ~ 0.2iRabi case. 



Ultra-short adiabatic production of squeezed states 

As described in Ref. [l2| for a particle in a harmonic 
oscillator and ultrafast adiabatic evolution, the method 
presented requires a certain period of the evolution where 
the trapping potential is actually reversed into an ex- 
pelling potential (an inverted parabola). Through the 
previously explained mapping, this corresponds in our 
many-body BH problem to a period of time where the 
atom-atom interaction is switched from repulsive to at- 
tractive, A < 0. For attractive interactions the poten- 
tial V{z) becomes a double- well potential in Fock space 
and the parabolic approximation of Eq. ([3]) does only 
hold for central regions in z [l7| . Thus, one can expect 
that the present method should still work only when the 
wavefunction tp^Zjt) does not spread in z significantly 
during the time where the atom-atom interactions are 
attractive. As seen in Fig. [3] this is indeed the case. We 
consider tp = O.OStRabi, A(G) = and A{tp) = 500. 
Then the required A{t) (upper panel) becomes negative 
for 0.025 < t/*Rabi ^ 0.04. The method is seen to work 
accurately, producing a fidelity ~ 1 at the final time. 
The intermediate evolution is however highly diabatic. 
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FIG. 4. (color online) Coherent spin squeezing attainable 
using a simple "Bang-bang" method as described in the text 
(lower panel). The fidelity between the evolved state and its 
corresponding instantaneous evolution is given in the upper 
panel. A(0) = in all cases, while A(fbang) ~ 5 (dashed), 20 
(dotted), and 200 (solid). N = 100. 



with close to zero fidelity between the evolved state and 
its corresponding adiabatic one, during the A(t) < 
evolution. The correspondence between the many-body 
BH and the single particle problem is made through the 
model of Eq. ([2]), which, as in all the other cases pre- 
sented agrees well with the BH results. 

Another set of protocols known to produce fast- 
adiabatic evolution in harmonic oscillator Hamiltonians 
arc the so-called "bang-bang" methods. The transla- 
tion of one of the simplest versions of these methods 
to the present problem requires the following steps as- 
suming that we start as before from an initial ground 
state corresponding to A(0) = Aq, and want to end up 
in the ground state corresponding to A^?: a) at i = 
set A = Ab = uJoutp/'i: — 1 (where ujq = ^/l + Aq and 
Lup = \/l + Ai?), b) let the system evolve at fixed Ab 
until the time t = Ib = ^Rabi/ {"^y/^^o^F), and c) change 
Ab to Ai?. Note that the final time is not arbitrary, as 
before, but is fixed by the initial and final values of A. 
This makes it in practice more unstable with respect to 
small errors [3] ■ Under ideal conditions (perfect timing) 
the method is seen to work well producing fairly low val- 
ues of ^1 at times of the order of ^Rabi/lO for iV = 100 
atoms with a fidelity of ~ 1 at the final time, see Fig. 21 

In this letter we have presented protocols for fast gen- 
eration of very coherent-spin-squeezed states in bosonic 
Josephson junctions. The attained squeezing is the one 
corresponding to the ordinary adiabatic evolution in the 
case of repulsive atom-atom interactions, but requires 
much shorter preparation times. Ordinary adiabatic 
squeezing is known to improve as N is increased as 
^1 oc [25j, faster than the single-axis-twisting Q 

^1 oc iV~^/^ or the recently proposed production of 



Thus, practical methods of fast-adiabatic 
driving have important advantages for any future im- 
plementation of BJJs where they are used as resources 
of entanglement /spin squeezing. The present procedures 
requires a good control of the time variation of the atom- 
atom scattering length during the desired period, a pos- 
sibility now at hand in current experimental setups for 
internal BJJ's. The methods have been obtained by ex- 
tending recently developed protocols for fast-adiabatic 
evolution in the case of a single particle evolving in a time 
dependent Harmonic oscillator potential, to the Bose- 
Hubbard Hamiltonian. If confirmed experimentally, the 
present proposal represents a useful step forward in the 
fast-adiabatic preparation of many-body entangled quan- 
tum resources. 
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